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(57) ABSTRACT 

A method for multi -objective fault accommodation using 
predictive modeling is disclosed. The method includes using 
a simulated machine that simulates a faulted actual machine, 
and using a simulated controller that simulates an actual 
controller. A multi -objective optimization process is per- 
formed, based on specified control settings for the simulated 
controller and specified operational scenarios for the simu- 
lated machine controlled by the simulated controller, to gen- 
erate a Pareto frontier-based solution space relating perfor- 
mance of the simulated machine to settings of the simulated 
controller, including adjustment to the operational scenarios 
to represent a fault condition of the simulated machine. Con- 
trol settings of the actual controller are adjusted, represented 
by the simulated controller, for controlling the actual 
machine, represented by the simulated machine, in response 
to a fault condition of the actual machine, based on the Pareto 
frontier-based solution space, to maximize desirable opera- 
tional conditions and minimize undesirable operational con- 
ditions while operating the actual machine in a region of the 
solution space defined by the Pareto frontier. 

17 Claims, 8 Drawing Sheets 
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METHOD AND SYSTEM FOR FAULT 
ACCOMMODATION OF MACHINES 

FEDERAL RESEARCH STATEMENT 

5 

This invention was made with Government support under 
contract NAS3-01 135-Task#3 awarded by NASA. The Gov- 
ernment has certain rights in this invention. 

BACKGROUND OF THE INVENTION 10 

The present disclosure relates generally to fault accommo- 
dation in complex engineered systems using predictive mod- 
eling and optimization. 

Fault accommodation is becoming an increasingly impor- 
tant subject for complex, fault-tolerant engineered systems, 
and in particular aircraft engine systems. Faults, if not suit- 
ably responded to in a timely manner, may lead to undesirable 
scenarios, negatively affecting operational safety. In aircraft 2 o 
engines there are a large number of fault root causes. Some of 
these faults in the high-pressure compressor (HPC) and high- 
pressure turbine (HPT) modules have symptoms that closely 
resemble highly deteriorated component states. Such faults 
may be accommodated by a structured manipulation of the 25 
engine control systems. 

However, due to the highly nonlinear nature of the engine 
controller and the fact that it is implemented as a large col- 
lection of computer modules (typically over 1 00) that employ 
a variety of one- and two-input tables, switching variables, 30 
logical elements, limiters, and priority-select logic, to name a 
few, the control design space is high-dimensional, highly 
nonlinear, multimodal, and discontinuous. To find an optimal 
accommodation, it is very important, yet non-trivial, to define 
the performance metric in a flexible and non-analytical man- 35 
ner. This is necessary in order to properly account for such 
diverse requirements as maintaining stall margins above cer- 
tain limits, minimizing both peak temperatures and the time 
spent above a certain temperature, and obtaining short rise 
times in response to changes in demand values. Furthermore, 40 
the changes must be accomplished over a wide range of flight 
conditions and disturbance inputs. 

Only a very small portion of an overall engine control 
system is designed to operate in a linear fashion, and even 
then, the controller gains are often scheduled as functions of 45 
the operating conditions (altitude, Mach number, and ambi- 
ent temperature deviation from standard day, for example). 
Although much is known about the behavior and design of 
linear control systems, this information is not relevant to the 
problems under consideration here. Rather, one must be pre- 50 
pared to work in the nonlinear domain, where theories and 
analytical results are much more scarce than for the linear 
domain. Also, the literature on nonlinear control systems, of 
necessity, tends to deal with specific situations, such as the 
area of integrator-windup protection (IWP). 55 

It is not to be expected that conventional optimization 
methods and those that depend on gradient evaluations should 
work, and in view of the existing art, it would be beneficial to 
provide an application of evolutionary algorithms to aircraft 
engine control systems design, where the controls design and 60 
optimization is performed using a lull-order engine model 
and full control systems structures that do not oversimplify 
the inherent complexities in these highly complex nonlinear 
dynamic systems. Accordingly, there is a need in the art for a 
non-conventional optimization method to provide real-time 65 
adj ustments to the controller to recover from some or all of the 
engine faults. 
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BRIEF DESCRIPTION OF THE INVENTION 

An embodiment of the invention includes a method for 
multi-objective fault accommodation using predictive mod- 
eling. The method assumes the presence of a diagnostic rea- 
soner that indicates the presence of a fault. The method 
includes using a simulated machine that simulates a faulted 
actual machine, and using a simulated controller that simu- 
lates an actual controller, the simulated machine being con- 
trolled by the simulated controller, and the actual machine 
being controlled by the actual controller. A set of multi - 
objective optimization processes are performed which 
include 1) an off-board (off-line) process for pre-generating 
safety-restoring optimal conditions for quick lookup retrieval 
when a fault is detected, and 2) either an on-board (on-line) or 
an off-board (off-line) process that calculates optimal solu- 
tions with an expanded set of objectives. Each process is 
based on specified control settings for the simulated control- 
ler and specified operational scenarios for the simulated 
machine controlled by the simulated controller, to generate a 
Pareto frontier-based solution space relating performance of 
the simulated machine to settings of the simulated controller, 
including adjustment to the operational scenarios to represent 
a fault condition of the simulated machine. Control settings of 
the actual controller are adjusted, represented by the simu- 
lated controller, for controlling the actual machine, repre- 
sented by the simulated machine, in response to a fault con- 
dition of the actual machine, based on the Pareto frontier- 
based solution space, to maximize desirable operational 
conditions and minimize undesirable operational conditions 
while operating the actual machine in a region of the solution 
space defined by the Pareto frontier. 

An embodiment of the invention includes a system for 
multi-objective fault accommodation using predictive mod- 
eling. The system includes a simulated machine that simu- 
lates a faulted actual machine, a simulated controller that 
simulates an actual controller, the simulated machine being 
controlled by the simulated controller, and the actual machine 
being controlled by the actual controller, a processor, and an 
adjuster portion. The processor performs a multi-objective 
process, based on specified control settings for the simulated 
controller and specified operational scenarios for the simu- 
lated machine controlled by the simulated controller, to gen- 
erate a Pareto frontier-based solution space relating perfor- 
mance of the simulated machine to settings of the simulated 
controller, including adjustment to the operational scenarios 
to represent a fault condition of the simulated machine. The 
adjuster portion adjusts control settings of the actual control- 
ler, represented by the simulated controller, for controlling 
the actual machine, represented by the simulated machine, in 
response to a fault condition of the actual machine, based on 
the Pareto frontier-based solution space, to maximize desir- 
able operational conditions and minimize undesirable opera- 
tional conditions while operating the actual machine in a 
region of the solution space defined by the Pareto frontier. 

An embodiment of the invention includes a computer read- 
able medium for multi -objective fault accommodation using 
predictive modeling, the computer readable medium includes 
computer executable instructions for facilitating an embodi- 
ment of the aforementioned method. 

These and other advantages and features will be more 
readily understood from the following detailed description of 
preferred embodiments of the invention that is provided in 
connection with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Referring to the exemplary drawings wherein like elements 
are numbered alike in the accompanying Figures: 
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FIG. 1 depicts in block diagram form a computational 
model having sub-blocks in accordance with embodiments of 
the invention; 

FIG. 2 depicts in schematic form a proportional path in a 
VSV actuator control system in accordance with embodi- 5 
ments of the invention; 

FIG. 3 depicts in schematic form a proportional path in a 
FMV actuator control system in accordance with embodi- 
ments of the invention; 

FIG. 4 depicts in block diagram form interconnected con- io 
trol subsystems in accordance with embodiments of the 
invention; 

FIG. 5 depicts in graph form results of accommodation for 
several performance parameters over time in accordance with 
embodiments of the invention; 15 

FIG. 6 depicts in graph form an exemplary Pareto frontier- 
based solution space having a Pareto frontier in accordance 
with embodiments of the invention; 

FIG. 7 depicts in block diagram form a system for practic- 
ing methods in accordance with embodiments of the inven- 20 
tion; and 

FIG. 8 depicts a flow diagram for a system in accordance 
with embodiments of the invention. 

DETAILED DESCRIPTION OF THE INVENTION 25 

Embodiments of the invention relate to fault accommoda- 
tion in complex engineered systems such as, but not limited 
to, aircraft engines, gas turbines, mechanical systems, chemi- 
cal processing systems, and electromechanical systems, each 30 
with dedicated control systems. 

An aircraft engine is composed of interconnected mechani- 
cal parts and a computer control system that coordinates their 
operation. Faults are typically distinguished from gradual 
deterioration. They are more severe in magnitude and occur at 35 
a timescale with a progression much faster (by an order of 
magnitude) than ordinary wear. Ordinarily, faults cannot be 
predicted because they are due to factors that cannot directly 
be observed. An example is inhomogenieties in material that 
over time may lead to cracks in a blade. The method described 40 
below assumes that a diagnostic reasoner indicates the pres- 
ence of a fault. 

Computer simulation models for the engine and its accom- 
panying control systems are described and employed herein. 

In order to be able to perform fault accommodation on the 45 
actual engine on wing, first it is necessary to evaluate and 
identify what compensation strategies will work in simula- 
tion. By having a good sense of what compensatory method 
will work, as identified through a simulation-based evalua- 
tion, it is possible to deploy that same strategy on the engine 50 
control system while the engine is on wing. As such, it is 
highly desirable that the engine and control system simula- 
tion models be highly reliable and have a high response 
fidelity. 

Fault accommodation is performed at the engine control 55 
systems level by adjusting various controller characteristics. 
This accommodation is done at the controller level since 
mechanical systems are not adaptable to change in the man- 
ner of computer control software settings. Accommodation at 
a mechanical level would require at least partial disassembly 60 
of the engine (either by humans or robots) and maintenance 
(either repair or replacement), which is time consuming and 
expensive. By making changes to the engine control system 
settings to achieve the same or better result, undesirable main- 
tenance time and expense may be avoided. Of course, for 65 
serious problems, or during scheduled maintenance intervals, 
mechanical disassembly and fix would be performed. 
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In order to perform fault accommodation, as disclosed 
herein, the simulation is operated at the level of engine-plus - 
controller. Therefore, and for all practical purposes, the simu- 
lation is representative of the actual engine with its controller. 
An optimizer, specifically a multi-objective optimizer, is 
interfaced to the simulation models to identify those control 
settings that result in the most desirable compensated system 
behavior. A multi-objective optimizer is utilized to simulta- 
neously consider multiple desirable compensated system 
behavioral measures, and help identify the best tradeoffs in 
this system behavioral space. 

In an aircraft engine, multiple measures of dissatisfaction 
of system behavior may be present, which results in a goal to 
simultaneously minimize multiple dissatisfaction measures. 
One example of a dissatisfaction measure is the distance away 
from ideal new engine performance measures, such as stall 
margin, thrust, and peak exhaust gas temperature, for 
example. To minimize these dissatisfaction measures, a 
Pareto frontier is identified, with each point on the Pareto 
frontier tying back to a particular control adaptation strategy. 
Subsequent to establishing a Pareto frontier-based solution 
space, a particular point on the Pareto can be down- selected 
using a variety of methods, where a corresponding control 
adaptation strategy can be selected and then deployed to the 
actual engine’s controller to obtain the desired compensated 
behavior. 

In view of the complexity of an aircraft engine and aircraft 
engine controllers, it is understood that one skilled in the art 
is a person having knowledge of prior art aircraft engine and 
aircraft engine controller architectures, and therefore details 
of such architectures are not presented herein. 

An exemplary evolutionary search algorithm is disclosed 
herein having features and characteristics that make it par- 
ticularly well suited to the problem of HPC and HPT fault 
accommodation. More background information is presented 
on multi -objective evolutionary algorithms and their applica- 
tion to aircraft engine control systems design. 

An embodiment of the invention discloses a method to 
automatically adjust controllers for aircraft engines such that 
HPC and HPT faults might be temporarily accommodated 
in-flight. How the optimization of different elements within 
an overall controller can be addressed is shown in an efficient 
fashion. These elements include local actuator gains, control 
modifiers, and control schedules. An evolutionary algorithm, 
which will be discussed in more detail below, is utilized to 
realize multi-objective optimization on a local as well as a 
global level, depending on the optimization task at hand. 
Fitness functions comprise performance metrics that incor- 
porate stall margins, exhaust gas temperature, fan-speed 
tracking error, and local tracking errors. 

To illustrate the challenges involved with making changes 
to the controller, the controller design problem is first illumi- 
nated. Typically, aircraft engine controller design is an itera- 
tive process. Initially, a linear engine model is built by extract- 
ing partial derivatives from models based on first-principles 
(Close et al., 2001). Then, local controllers are designed and 
optimized using first-principles as well as derivatives from 
previous engine designs. Next, schedules are designed using 
performance requirements. Finally, the control logic is estab- 
lished which integrates the individual components and tales 
overall stability and performance requirements into account. 

The performance of the overall control system is tested on 
increasingly more complex systems starting with the local 
model, the bare component level model (CLM), the CLM 
with the full controller integrated, a dry rig test, wet rig test, 
test cell runs, and test flight. Each test cycle might necessitate 
a revision of some controller components with renewed vali- 
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dation and verification. There are a number of different cat- 
egories which are affected by the design and which could also 
be considered for accommodation: 

local actuator gains, either constant or scheduled, 
logic thresholds, 

adders and multipliers for gains and schedules, 
schedule entries, and 
control logic structure 

For embodiments of the invention disclosed herein, opti- 
mization is performed for select control variables in the first 
four categories. However, design changes in the control logic 
structure are contemplated. 

Optimization of the controller for nominal engines is bro- 
ken up into several complementary subtasks. These subtasks 
include: (i) optimization of the actuator gains, (ii) optimiza- 
tion of the control modifiers (adjustables), and (iii) design and 
optimization of the control schedules. This task decomposi- 
tion is a consequence of the fact that local gain modifications 
often do not result in any significant variation at the global 
performance level. In addition, the potential for crosstalk, that 
is, the difficulty to track correlations of several simulta- 
neously manipulated variables on the overall controller, sup- 
ports the strategy of dividing the optimization endeavor into 
smaller optimization tasks. Depending on the impact the par- 
ticular control variable under consideration has on the overall 
and local performance criteria, we maximize the observabil- 
ity from an optimization standpoint. This means that for some 
control variables only local performance criteria (local track- 
ing errors, for example) are considered while other control 
variables are considered from a global level (stall margins 
(SM), exhaust gas temperature (EGT), and fan-speed track- 
ing error (nj, for example). 

FIG. 1 gives an overview of this strategy in the form of a 
model 100 having several sub-blocks. Information from sen- 
sors 130 along with desired performance input 105 are 
entered into a global and local performance metrics block 
110, the output of which is acted upon by an optimization 
module 115. 

Optimization can take advantage of a computational model 
simulator 100, herein referred to as an FSIM, which can 
simulate the dynamic behavior of a production aircraft engine 
and its controller with a high degree of fidelity. The simula- 
tion modules comprise the CLM 127 and an emulation of the 
Full Authority Digital Electronic Control (FADEC) 125. A 
user may specify control settings and flight scenarios 120, and 
execute FSIM to obtain the engine response given a high- 
level (pilot) command such as demanded fan speed, which is 
a good measure of thrust. In an embodiment, the pilot’ s thrust 
specification is via the position of the throttle resolver angle 
(TRA) measured in degrees, where a lower TRA setting cor- 
responds to a lower thrust demand. Sensors 130 measure 
various engine response characteristics due to an input thrust 
demand specification. In an embodiment, sensors 130 in com- 
bination with information derived by global and local perfor- 
mance metrics block 110 provide the functionality of a rea- 
soner for indicating the presence of a fault. 

Prior to the invention disclosed herein, and for a product 
development project looking into automated controller 
adjustments for deteriorated engines, a combination of 
domain knowledge and several trial runs of FSIM have been 
utilized to isolate and identify actuator gain candidates to be 
optimized. Based on the principal consideration of the impact 
a particular actuator has on engine performance, and fault 
accommodation, the Fuel Metering Valve (FMV) propor- 
tional gain (FMV-Kp), and the Variable Stator Vane (VSV) 
proportional gain (VSV-Kp) were selected as parameters to 
be optimized. However, it is appreciated that other actuators 
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can be chosen for fault accommodation as well. While VSV- 
Kp is a constant, FMV-Kp is a function of the FMV actuator 
position and the corrected core speed, where a higher core 
speed results in a higher gain. A simplified view of the pro- 
5 portional path 135 in the VSV actuator control system is 
shown in FIG. 2, and a simplified view of the proportional 
path 140 in the FMV actuator control system is shown in FIG. 
3. 

In the actuator loops (proportional paths 135, 140) of 
10 FIGS. 2 and 3, d is the demanded actuator position, and p is 
the achieved actuator position. The optimization problem 
then is to select that gain value that minimizes the time inte- 
gral of the square of the position tracking error e of the 
actuator loop, which is represented by the following equation. 


min./ = J e 2 dX 

20 

The actuator position demand signal d is computed by the 
FSIM simulation modules through a complex transformation, 
known to one skilled in the art, of the thrust demand profile 
(with respect to time) specified by the pilot. 

25 Referring now to FIG. 4, the engine control system is large 
and complex, but known to one skilled in the art, with numer- 
ous interconnected subsystems such as the FMV 145 and 
VSV 150 subsystems. Each of these control systems is pro- 
vided with a suite of modifiers 155, alternatively also referred 
30 to as adjustables, that consists of adders and multipliers for 
gains, adders and multiplier for schedules, and logic thresh- 
olds. While only the FMV 145 is illustrated in expanded block 
diagram form in FIG. 4, it will be appreciated that the VSV 
150 and other control subsystems will have a similar 
35 expanded block diagram form. 

In an embodiment, the FMV control system has a set of 53 
modifiers, while the VSV control system has a set of 20 
modifiers. Each modifier is adjustable and is a bounded real 
number, and the bounds are specified in FSIM. The optimi- 
40 zation problem is to identify a set of modifiers such that global 
performance criteria are optimized. Ideally, to attain robust- 
ness in solution quality, the performance metric should be 
cumulatively considered over a number of flight conditions 
such as: 

45 altitude, 

Mach number, and 

ambient temperature deviation from standard day and 
engine configurations such as: 
customer bleed, 

50 horsepower extraction, 
deterioration, and 
component tolerances. 

To adequately evaluate performance, a metric is applied 
that allows the quantification of all global performance 
55 requirements. Such an ideal global performance metric may 
include one or more of the following relative measures: 
booster stall margin vs booster inlet flow, 
compressor stall margin vs compressor inlet flow, 

VSV demanded position vs corrected core speed, 

60 VBV demanded position vs corrected fan speed, 
corrected Phi vs corrected core speed, 
combustor fuel/air ratio vs severity parameter, 
high-pressure turbine inlet temperature vs corrected core 
speed, and 

65 exhaust gas temperature vs corrected fan speed and time. 
As used herein, Phi is the ratio of the fuel flow and the 
corrected combustor static pressure. 
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To reduce the complexity of the global performance met- 
ric, we focus on typical input variations and study the most 
important parameters for aircraft engine control systems vali- 
dation. In particular, we drive to meet all stall margin limits, 
good tracking of a fan- speed demand profile, and reduction in 
the peak exhaust gas temperature. Towards this end, let E be 
the exceedance profile comprising the EGT exceedance 
E EGr> the fan stall margin exceedance E w , the booster stall 
margin exceedance Esm 2 > and the compressor stall margin 
exceedance Ew 25 3 and EGT be the exhaust gas temperature 
profile, EGT min the table cruise temperature, e the exponent 
by which the distance to the limit is penalized, a a weight for 
the temperature component, b a weight for the fan-speed 
tracking error component, n l the fan-speed profile, n x the 
fan-speed demand profile, t the time. Then the multi-objective 
optimization problem has the vectorial form min J, where 


J = |f; J max(0, a ■ ( EGT - EGT min ) e ) dX\ bj \n { - n l<Jmd \ dX 
and 


ance E SMl2 , the booster stall margin exceedance E SMi , and the 
compressor stall margin exceedance E SM ^, n x the fan-speed 
profile, n^the fan-speed demand profile described as a step 
function with the step coinciding with the burst phase of the 
5 Bodie, t the time. Then the multi -objective optimization prob- 
lem has the vectorial form min J, where: 


10 


f l«i 

t 



and 




E = E S M n + i 

£571*2 + £571*25 


Esm 12 = " 

f 00 

if SM l2 <SM l2min 

15 

) 



lo 

otherwise 


( 

00 

if SM 2 <SM 2min 


ii 

0 

otherwise 

20 

£571*25 — 

j 00 

lo 

if SM 25 <SM 25min 
otherwise 


E = Eegt + E S m 12 + Esm 2 + Esm 25 


Eegt = j 

0 

00 

if EGT < EGT max 
otherwise 

25 

Esm 12 = • 

(7 

if SM l2 < SM i2 min 
otherwise 


II 

§ 

00 

0 

if SM 2 < SM 2min 
otherwise 

30 

II 

(7 

if SM 25 < SM 2 5min 
otherwise 



The control logic in a typical aircraft engine controller 
utilizes a suite of schedules that are functions of one or two 
input variables. Schedules are typically but not necessarily 
implemented as lookup tables and the output values are com- 
puted via linear interpolation among the closest neighbors. 40 
Schedule surfaces (output maps) represent nonlinear trans- 
formations of the inputs to the output and are important com- 
ponents of an aircraft engine’s control logic. 

Based on a combination of domain knowledge, simulation, 
and knowledge elicitation from domain experts, a particular 45 
control schedule, here denoted the F136 schedule (FI 36 is a 
control schedule within an aircraft engine available from 
General Electric Company) in the FMV module was selected 
as a candidate for use herein for evolutionary optimization. 
This schedule outputs a rate-gain reduction given the ambient 50 
pressure (a physical function of altitude) and compressor 
speed, and is active during the burst phase for a specific 
aircraft maneuver called a Bodie, wherein at cruise the pilot 
cuts thrust for a short time period and increases thrust through 
a burst before the engine temperatures have achieved steady 55 
state at the reduced power level. In the absence of a rate-gain 
schedule, the fan-speed response during the burst phase is 
extremely sluggish, which is highly undesirable. What is 
desirable however, is a rapid return to the original fan speed. 

In an embodiment, an optimization problem then is to 60 
identify a set of F 136 schedule entries such that fan accelera- 
tion is maximized during the burst phase of the Bodie, subject 
to maintaining all stall margins above acceptable limits. 
Simulation-based evaluation reveals that during this maneu- 
ver the EGT is always well within limits and is therefore not 65 
included in the global performance metric. Let E be the 
exceedance profile comprising the fan stall margin exceed- 


An evolutionary optimization of a schedule that is the 
function of two input variables corresponds to a systematic 
and joint manipulation of the table entries. An important 
aspect in the optimization of these control surfaces is the 
smoothness of these derived surfaces. Unless smoothness is 
explicitly included as a design requirement, an evolutionary 
optimization can result in noisy, albeit optimal, schedules. To 
facilitate smoothness in derived schedule surfaces, the entries 
in each test surface T are filtered using a specialized bi- 
directional filtering algorithm that is applied to each derived 
row T (i, j), and is shown below. 


£*(*»= o i 

TJh j)=aT A (iJ+l)+(l j) 2 

T s (i,-l)=aT A (i,0) 3 

T& j)=aT s (i, j-l)+(l-a)T(i, j ) 4 


In the algorithm presented above, a is a smoothing factor. 
The objective of the algorithm is to first identify a reliable 
starting value T s (h -1) (line 3) for each row i following the 
procedure outlined in lines 1, 2, and 3. This is an important 
step, since a quality outcome is dependent upon selection of a 
reliable starting point. Next, the smoothed values T s (i, j) are 
computed using the procedure outlined in line 4. 

In an embodiment, and for simulation and optimization 
purposes, the effect of a faulty engine can be modeled by 
decreasing the efficiencies and flow scalars of the rotating 
components (as an aside, it is noted that the symptoms of 
deteriorated engines and engines with some non-catastrophic 
HPC/HPT faults are similar except that the rate of change of 
progression of symptoms at fault initiation is often times an 
order of magnitude higher than the rate of change of deterio- 
ration related symptoms). The major modules affected for a 
commercial high-bypass engine are the fan, booster, com- 
pressor, high-pressure turbine, and low-pressure turbine. 
Relative adjustments to these variables are shown in Table 1 , 
which also shows for comparison adjustments that would be 
made for small and large HPT and HPC faults for an exem- 
plary commercial, high-bypass, twin-spool, turbofan engine, 
respectively. 
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TABLE 5 


Typical Fault and Deterioration Adjustments 
for an exemplary engine simulator 



HPC 

HPC 

HPT 

HPT 


Efficiency 

Flow 

Efficiency 

Flow 


Scalar 

Scalar 

Scalar 

Scalar 

Small HPT Fault 



-1.50% 

0.15% 

Large HPT Fault 
Small HPC Fault 

-1.50% 

-1.50% 

-6% 

0.60% 

Large HPC Fault 

- 6 % 

-6% 




For deterioration, other adjustments are made to the fan 
module, booster module, and LPT module to capture the 
system wide wear. Typical adjustments are shown in Table 2 
for a different engine type. 


TABLE 2 


Typical Deterioration Adjustments 


Efficiency 


Component 

Scalar 

Flow Scalar 

Fan 

-0.015 

-0.5% 

Booster 

-0.001 

-0.6% 

LPT 

-0.011 

+0.4% 


Since the engine controller is designed to follow the pilot’s 
demanded fan speed to closely meet specified thrust require- 
ments, the changes in a deteriorated (and otherwise not 
faulted) engine result in higher fuel consumption and higher 
temperatures of the high- and low-pressure turbine blades. 
Table 3 gives an example of changes from a nominal engine to 
a large HPT and HPC fault . For comparison, we list again 
also signatures for a 50% deteriorated engine. Tables 3-5 
show small, and large HPC and HPT faults in comparison to 
50% deterioration. 


TABLE 3 


Typical Fault and Deterioration Effects 


% Delta 

% Delta 

% Delta (50% 

Sensor 

(Large HPT Fault) 

(Large HPC Fault) 

Deterioration) 

T12 

-0.37% 

-0.37% 

0.3% 

XN1 

0.00% 

-1.7% 

-0.50% 

XN2 

-0.89% 

-1.3% 

-1.25% 

T25 

2.94% 

4.2% 

0.25% 

P25 

2.35% 

1.7% 

0.19% 

T3C 

-2.60% 

3.2% 

-1.08% 

PS3 

-1.00% 

-5.2% 

-1.92% 

EGT 

1.32% 

5.5% 

1.23% 


TABLE 4 


Typical HPC Fault and Deterioration Effects 


% Delta 

% Delta 

% Delta (50% 

Sensor 

(small HPC Fault) 

(Large HPC Fault) 

Deterioration) 

T12 

-0.37% 

-0.37% 

0.3% 

XN1 

-0.30% 

-1.7% 

-0.50% 

XN2 

-0.3% 

-1.3% 

-1.25% 

T25 

1.0% 

4.2% 

0.25% 

P25 

0.5% 

1.7% 

0.19% 

T3C 

0.8% 

3.2% 

-1.08% 

PS3 

-1.3% 

-5.2% 

-1.92% 

EGT 

1.3% 

5.5% 

1.23% 


Typical HPT Fault and Deterioration Effects 


Sensor 

% Delta 

(small HPT Fault) 

% Delta 

(Large HPT Fault) 

% Delta (50% 
Deterioration) 

T12 

-0.37% 

-0.37% 

0.3% 

XN1 

-0.50% 

-1.8% 

-0.50% 

XN2 

-1.7% 

-4.4% 

-1.25% 

T25 

1.1% 

4.1% 

0.25% 

P25 

0.5% 

1.7% 

0.19% 

T3C 

-1.6% 

-3.9% 

-1.08% 

PS3 

-2.3% 

-8.3% 

-1.92% 

EGT 

1.7% 

6.9% 

1.23% 


15 It is seen that the HPC and HPT faults and deterioration are 
coupled to quite a degree. Not shown here is a common 
reduction in stall margins and an increase in thrust and fuel 
consumption. 

In an embodiment, changing controller behavior is further 
20 developed by selecting the modifiers of the FMV control 
logic and the modifiers of the VSV control logic as potential 
candidates for optimization to recover performance from a 
faulty engine. These modifiers are same ones described above 
in connection with FIG. 4. An evolutionary search algorithm 
25 is employed to find a set of FMV and VS V modifiers that meet 
all performance criteria. Again, the performance criterion 
used is as described in above in connection with FIG. 4. 

Discussion now turns to initial observations followed by 
systematic experiments that led to an embodiment of the 
30 simulation-based optimization of actuator gains, controls 
modifiers, schedules, and control modifiers for faulted 
engines. 

In selecting candidates to which the generic search algo- 
rithm was to be applied, several considerations were taken 
35 into account. For one, the search was started in an area of 
constrained complexity. The first subject was the code that 
handles the adjustments to the variable stator vane (VSV) 
positions. This operation can be divided into two parts: (i) the 
determination of the demanded VSV position, based on cur- 
40 rent flight conditions; and, (ii) the gains and other parameters 
of the VSV actuator loop itself. 

Upon looking at Beacon diagrams (Beacon is a computer 
program available from Applied Dynamics, Inc. that repre- 
sents an aircraft engine controller in terms of block diagrams 
45 and computational flow diagrams from which the actual com- 
puter code is generated) that describe the code used in the 
ECU (electronic control unit), it was concluded that the gains 
of the VSV actuator loop made a good starting starting point 
to demonstrate this approach. Specifically, experiments were 
50 carried out with the proportional gain of that loop, since the 
code had an adjustable adder defined which could be set to 
any desired value at the start of a rum, without having to 
rebuild the FSIM code. 

Several FSIM runs were made with preset values of VSV- 
55 Kp for a combination of a burst (increase of TRA from 36 to 
78 degrees), a constant TRA for 25 seconds, followed by a 
chop (decrease in TRA from 78 to 36 degrees). For a given 
FSIM run, the value of the proportional-gain adder was varied 
such that the actual proportional gain was a constant between 
60 the limits of 78 and 778. An exemplary design value is 578, 
and it is a constant (not scheduled according to core speed or 
any other variable). In terms of the major engine variables, 
such as fiiel flow, fan speed, VSV angle, and exhaust gas 
temperature (EGT), for example, there were virtually no dif- 
65 ferences over the 10-fold range of the proportional gain. The 
only variable that was affected by the change in proportional - 
gain value was the error in the VSV actuator loop. This 
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observed behavior suggests that the proportional gain could 
be selected so as to minimize the square of the actuator error. 

Optimization of the VSV-Kp follows the procedure out- 
lined above in connection with FIGS. 2-3. Engine operation 
was simulated subject to changes in throttle position while 
cruising at 35,000 ft., Mach 0.8, and standard-day tempera- 
ture. 

A burst and Bodie were used to excite the overall system, 
and gain optimization was performed independently for each 
of these excitation profiles. While no attempt was made to 
determine optimized values for other gains in the loop, such 
as the integral gain, or other parameter values, such as a 
lead-time constant, it is contemplated that it would certainly 
be possible to do so. 

Next, and with respect to runs with constant values of 
FMV-Kp, an embodiment of the actuator design method was 
applied to the fuel metering valve (FMV) actuator loop pro- 
portional gain. Here, the actuator controller is slightly more 
complicated than the VS V one in that the proportional gain is 
a function of both the FMV actuator position and the core 
speed, where a higher core speed results in a higher gain. As 
in the VSV actuator study, it was found that changes in the 
FMV proportional gain had virtually no effect on the engine 
variables (fan speed, fuel flow, and temperatures, for 
example). The scheduling of the gain was removed by modi- 
fying the appropriate schedule table (FI 32) and the burst and 
chop simulations were run over a wide range of constant 
gains. As before, only the error in the actuator loop was 
affected. 

Optimization of the FMV-Kp follows the procedure out- 
lined above in connection with the discussion relating to 
FIGS. 2 and 3. Engine operation was simulated subject to 
changes in throttle position while cruising at 35,000 ft., Mach 
0.8, and standard-day temperature. 

In an embodiment, a burst and Bodie were used to excite 
the overall system, and gain optimization was performed 
independently for each of these excitation profiles. In order to 
find important parameters whose values will affect global 
performance metrics (such as stall margins or exhaust gas 
temperature, for example), attention was focused on the 
extensive set of ECU modules that are used to produce the 
incremental changes in the demanded fuel flow. These incre- 
mental changes are continuously summed in order to produce 
the demand value for the FMV actuator. The behavior during 
a burst operation, in response to a sudden request from the 
pilot for an increase in power, was examined. During such a 
situation, several different regulators are active, depending on 
the various limits and constraints that must be satisfied in 
order to guarantee safe operation of the engine. By looking at 
which regulator was selected as time evolved following the 
burst command, the sequence of active regulators was deter- 
mined. 

Discussion is now directed to optimization of control sys- 
tem modifiers for a faulted engine. Engine operation was 
simulated over the range of faulted engines in an off-board 
optimization. The results were compiled into a lookup 
scheme, such as a lookup table for example, for quick or 
immediate retrieval. Operating points for which a direct 
match was not available were interpolated between existing 
solutions. The quick retrieval is best illustrated by an example 
case. At cruise for a nominal, healthy engine, a fault was 
inserted at 0.5 seconds and then 2 seconds later (reflecting the 
time a diagnostic tool might take to detect the fault), the 
adjustment parameters were changed to the accommodated 
variables determined by the off-board genetic algorithm. The 
runs continued for 20 seconds to allow transients to settle and 
achieve the steady- state accommodated measures for oper- 
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ability and performance. The modifiers considered here were 
VBV Adjustment Adder, VSV Adjustment Adder, Max Core 
Speed Adder, Horsepower Extraction, ACC Multiplier, ACC 
Adder, Throttle, Aux Bleed Switch, and Cowl Anti-Ice 
5 Switch. 

FIG. 5 [ 11 ] shows the temporal responses for this test case, 
illustrating Fan Stall Margin, Booster Stall Margin, Compres- 
sor Stall margin, and Net Thrust. Although the optimizer 
considered for this example was only a single joint objective 
to (a special and more limiting case of the multi-objective opti- 
mization scheme), all performance metrics recover quickly to 
safe conditions (illustrated by stabilized characteristic 
curves). Note that at this point (at the 20-second mark of 
simulation time) the fault itself still exists but its impact on the 
15 performance measures has been accommodated (illustrated 
by stabilized characteristic curves). 

The discussion above showed results for accommodation 
of the fault using precompiled solutions from off-board opti- 
mization. In a next step, it is contemplated that further objec- 
20 tives are being considered to arrive at a solution that provides 
optimality for a wider time window. Because the primary 
objectives relating to safe operation have been met, more time 
can be spent on secondary objectives, including the comple- 
tion of the mission (which could be the single flight or a set of 
25 flights), planning and scheduling (for example, airline sched- 
ules), resource availability to complete maintenance at a par- 
ticular site, risk of secondary damage, to name a few. The 
computation of solution for this expanded set of objectives 
can be accomplished either on-board or off-board. 

30 Also, in view of the many and varied solutions available, 
the results of an embodiment of the multi-objective optimi- 
zation disclosed herein may be expressed as a Pareto frontier, 
that is, a Pareto frontier-based solution space. In other words, 
there are a number of possible solutions that all meet the 
35 criteria that may result in different engine behavior. As such, 
it is then desirable to express a preference for a particular 
solution from the set of possible solutions. Ideally, the pref- 
erence could be cast within the objective function. However, 
the effects of the solutions are not always apparent until after 
40 the optimization results are reviewed. 

FIG. 6 illustrates an exemplary Pareto frontier-based solu- 
tion space 200 having a Pareto frontier 205 , which is plotted 
against two different engine characteristics, such as fuel con- 
sumption (Cf) and peak engine temperature (Tpk), for 
45 example. From FIG. 6, it will be appreciated that different 
solutions along the Pareto frontier 205 will result in different 
engine behavior, by virtue of the different fuel consumption 
and peak engine temperature characteristics. 

In summary, a framework to the accommodation of jet- 
50 engine controllers has been formulated by applying evolu- 
tionary search algorithms to actuators, multiplicative and 
additive adjustments, and table generation and/or modifica- 
tion. To that end, meaningful performance functions have 
been developed whose minimization produces controller 
55 parameters that result in desirable engine behavior. The meth- 
ods disclosed herein incorporate stall margins, peak EGT, and 
tracking of changing throttle positions. In addition, a smooth- 
ing function was employed that in effect penalizes discon- 
tinuous table solutions. Then, these techniques were used to 
60 successfully adjust up to 73 parameters at a time in the con- 
troller of a real commercial aircraft engine. In addition, the 
ability to tune the proportional gain of a regulator in the 
context of its operation in a nonlinear environment by mini- 
mizing the integral of the square of the actuator error was 
65 demonstrated. Moreover, evolutionary search algorithm 
methodology was employed to generate a 3-D table of the 
form z=f(x,y) to maintain rapid response to a demanded 
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power increase as a part of a Bodie maneuver. Finally, the 
methodology disclosed herein was applied to a faulted engine 
and showed that by adjustment of multiplicative and additive 
parameters, the engine’s stall-margin requirements can be 
met while also meeting other safety critical requirements. 

It is contemplated that the methodology disclosed herein 
will be useful for future work. One avenue is the integration of 
the optimization for design assistance during the various 
design and validation stages, cycle deck over CLM, FSIM, 
dry rig, wet rig, test cell, and flight test, for example. More- 
over, embodiment of the optimization approach disclosed 
herein may be extended to adapt engine performance based 
on in-service data, and to adapt engine performance as an 
engine deteriorates. This assistance could range from auto- 
mated validation of design choices to suggestion of param- 
eters as discussed and illustrated herein. Another avenue 
leads to scaling the optimization task. Of particular interest is 
ensuring cross -communication of individual results from 
components in a concurrent optimization scheme. Such a 
co-evolutionary optimization (Subbu and Sanderson, 2004) 
would allow concurrent module- specific exploration of the 
global design space, thus responding to the need of both 
domain- specific focus and adhering to global performance 
metrics, which could be accomplished via agent-based multi- 
objective optimization. Also of interest is the integration of 
external information such as expert knowledge, historical 
runs, and information from pilots during test flights, for 
example, which would need the development of an informa- 
tion aggregation component (Goebel et al., 2000, Goebel, 
2001) that can deal with the inherent uncertainties and the 
different format to more formally translate these observations 
into an objective function and performance metric. 

It is also contemplated that individualized optimizations of 
engines could be performed using modifiers by responding to 
specific engine characteristics (as opposed to model wide 
baselines), thus further improving performance. In addition, 
performance-enhancing optimization may be employed 
through the reduction of schedule size, thus reducing FADEC 
memory requirements and improving throughput, which 
could lead to a selection of optimal schedule size for a number 
of controllers such as the FMV and power management, 
which typically deal with large schedules. In addition, the 
overall FADEC architecture could be optimized by identify- 
ing obsolete elements (schedules, for example). Finally, the 
logic structure itself could be an opportunity for optimization, 
which could be accomplished through genetic programming 
or inductive learning such as Experience Based Learning. 

With respect to evolutionary algorithms (EAs), EAs 
include genetic algorithms (Goldberg, 1989, Holland, 1994), 
evolutionary programming (Fogel et al., 1966), evolution 
strategies (Back, 1996), and genetic programming (Koza, 
1992). The principles of these related techniques define a 
general paradigm that is based on a simulation of natural 
evolution. EAs perform their search by maintaining at any 
time t a population P(t)={P 1 (t), P 2 (t), . . . , P/>(t)} of individu- 
als. “Genetic” operators that model simplified rules of bio- 
logical evolution are applied to create the new and more 
superior population P(t+1). This process continues until a 
sufficiently good population is achieved, or some other ter- 
mination condition is satisfied. Each P z (t)€P(t) represents, via 
an internal data structure, a potential solution to the original 
problem. The choice of an appropriate data structure for rep- 
resenting solutions is very much an “art” than “science” due 
to the plurality of data structures suitable for a given problem. 
However, the choice of an appropriate representation is often 
an important step in a successful application of EAs, and 
effort is required to select a data structure that is compact, 
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minimally superfluous, and avoids creation of infeasible indi- 
viduals. For instance, if the problem domain requires finding 
an optimal real vector from the space defined by dissimilarly 
bounded real coordinates, it is more appropriate to choose as 
5 a representation a real-set-array (a real-set-array being an 
array of bounded sets of reals) instead of a representation 
capable of generating bit strings (a representation that gener- 
ates bit strings can create many infeasible individuals, and is 
certainly longer than a more compact sequence of reals). 

1 0 Closely linked to the choice of representation of solutions, 

is the choice of a fitness function J: P(t)^R, that assigns 
credit to candidate solutions. Individuals in a population are 
assigned fitness values according to some evaluation crite- 
15 rion. Fitness values measure how well individuals represent 
solutions to the problem. Highly fit individuals are more 
likely to create offspring by recombination or mutation opera- 
tions. Weak individuals are less likely to be picked for repro- 
duction, and so they eventually die out. A mutation operator 
20 introduces genetic variations in the population by randomly 
modifying some of the building blocks of individuals. Evo- 
lutionary algorithms are essentially parallel by design, and at 
each evolutionary step a breadth search of increasingly opti- 
mal sub-regions of the options space is performed. Evolution- 
25 ary search is a powerful technique of solving problems, and is 
applicable to a wide variety of practical problems that are 
nearly intractable with other conventional optimization tech- 
niques. Practical evolutionary search schemes do not guaran- 
tee convergence to the global optimum in a predetermined 
30 finite time, but they are often capable of finding very good and 
consistent approximate solutions. However, they are shown 
to asymptotically converge under mild conditions (Subbu and 
Sanderson, 2004). 

Most real-world optimization problems have several, often 
35 conflicting objectives. Therefore, the optimum for a multi- 
objective problem is typically not a single solution — it is a set 
of solutions that trade-off between objectives. The Italian 
economist Vilfredo Pareto first generally formulated this con- 
cept in 1 896, and it bears his name today. A solution is Pareto 
40 optimal if (for a maximization problem) no increase in any 
criterion can be made without a simultaneous decrease in any 
other criterion. The set of all Pareto optimal points is known 
as the Pareto frontier or alternatively as the efficient frontier. 

Pareto Frontier optimization techniques provide a frame- 
45 work for tradeoff analysis between, or among, desirable ele- 
ment attributes (e.g., where two opposing attributes for analy- 
sis may include turn rate versus range capabilities associated 
with an aircraft design, and the trade-off for an optimal turn 
rate (e.g., agility) may be the realization of diminished range 
50 capabilities). A Pareto Frontier may provide a graphical 
depiction of all the possible optimal outcomes or solutions. 
Evolutionary algorithms (EAs) may be employed for use in 
implementing multi-objective optimization functions. Multi - 
objective EAs involve searches for, and maintenance of, mul- 
55 tiple Pareto -optimal solutions during a given search which, in 
turn, allow the provision of an entire set of Pareto-optimal 
(Pareto Frontier) solutions via a single execution of the EA 
algorithm. 

A decision function may be applied to the Pareto Frontier 
60 for the decision-making selection. The decision function may 
be applied to the optimal sets of input-output vector tuples to 
reduce the number of input-output vector tuples in what may 
be referred to as a sub-frontier. One such decision function 
may be based on the application of costs or weights to objec- 
65 fives, whereby a subset of Pareto optimal solutions closest to 
an objectives weighting may be identified. Additional deci- 
sion functions such as one that is capable of selecting one of 
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the optimal input-output tuples that minimally perturbs the 
engine from its current state, may be applied. 

In view of the foregoing, it will be appreciated that embodi- 
ments of the invention include a method for performing 
multi-objective fault accommodation that uses a predictive 
model 100, based on specified control settings for a simulated 
controller 125 and specified operational scenarios for a simu- 
lated machine 127 controlled by the simulated controller, to 
generate a Pareto frontier-based solution space 200 relating 
performance of the simulated machine to settings of the simu- 
lated controller, including adjustment to the operational sce- 
narios to represent a non-catastrophic HPC (high pressure 
compressor) or HPT (high pressure turbine) fault condition of 
the simulated machine. With the model, control settings of an 
actual controller 125', represented and illustrated by the simu- 
lated controller 125, are adjusted for controlling an actual 
machine 127', represented and illustrated by the simulated 
machine 127, in response to a non-catastrophic HPC or HPT 
fault condition of the actual machine, based on the Pareto 
frontier-based solution space, to maximize desirable opera- 
tional conditions, such as low fuel consumption for example, 
and minimize undesirable operational conditions, such as 
high peak engine temperature for example, while operating 
the actual machine in a region of the solution space defined by 
the Pareto frontier 205. In an embodiment, the deteriorated 
condition of the actual machine is representative of normal 
wear of the actual machine. In an embodiment, the non- 
catastrophic HPC and HPT fault conditions of the actual 
machine are coupled. 

In an embodiment, the actual machine 127', and associated 
simulated machine 127, comprise interwoven subsystems 
320, 320' (see FIG. 7 for example), and the aforementioned 
multi-objective optimization process is performed in a feder- 
ated, or distributed, manner, thereby allowing the interwoven 
subsystems, in accordance with methodologies disclosed 
herein, to find their own respective solutions, while intercor- 
relations of the subsystems are controlled at a global, or 
system, level. In an embodiment, the aforementioned maxi- 
mizing and minimizing are performed off-line and the results 
stored onboard, or precompiled into lookup tables that are 
then stored onboard, such that they can be simply retrieved for 
the online case within a subsecond time frame sufficient to 
establish safety -restoring operating points for the actual 
machine 127', thereby accommodating the HPC or HPT fault 
condition of the actual machine 127'. In an embodiment, the 
aforementioned maximizing and minimizing are further per- 
formed in a time frame greater than a sub second subsequent 
to establishing the safety -restoring operating points to enable 
calculation of trade-offs on the Pareto frontier 205 to incor- 
porate contemporaneous demands on the actual machine 127' 
and external demands. In an embodiment, the contemporane- 
ous demands on the actual machine 127' include duration of 
operation of the actual machine, operational safety thresholds 
of the actual machine, deterioration of the actual machine, 
risk of secondary damage, logistic considerations for the 
actual machine, planning considerations for the actual 
machine, or any combination including at least one of the 
foregoing demands. 

In an embodiment, adjusting of the control settings include 
a bi-directional filtering algorithm, as set forth above, to 
facilitate smoothness in the derived schedule surfaces. 

In an embodiment, use of the predictive model, and adjust- 
ment of the control settings is performed on-board the actual 
machine, the aircraft for example, and the use of the predic- 
tive model and the adjustment of the control settings is per- 
formed at any time between consecutive operations of the 
actual machine, between consecutive flights of the aircraft for 
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example. In an embodiment, adjusting the control settings 
includes optimizing actuator gains, controls modifiers, and 
schedules of the CLM of the predictive model, and serve to 
reduce module temperature peaks, and to increase stall mar- 
5 gins of the aircraft engine. 

Also in view of the foregoing, and with reference now to 
FIG. 7, it will be appreciated that embodiments of the inven- 
tion also include a system 300 for multi -objective HPC or 
HPT fault accommodation. In an embodiment, the system 
10 300 includes a simulated machine 127 that simulates an HPC 
or HPT faulted actual machine 127', a simulated controller 
125 that simulates an actual controller 125', the simulated 
machine being controlled by the simulated controller, and the 
actual machine being controlled by the actual controller, a 
15 processor 305 that performs a multi-objective process, based 
on specified control settings for the simulated controller and 
specified operational scenarios for the simulated machine 
controlled by the simulated controller, to generate a Pareto 
frontier-based solution space 200 relating performance of the 
20 simulated machine to settings of the simulated controller, 
including adjustment to the operational scenarios to represent 
an HPC or HPT fault condition of the simulated machine, and 
an adjuster portion 310 that adjusts control settings of the 
actual controller 125', represented by the simulated controller 
25 125, for controlling the actual machine 127', represented by 
the simulated machine 127, in response to an HPC or HPT 
fault condition of the actual machine, based on the Pareto 
frontier-based solution space, to maximize desirable opera- 
tional conditions and minimize undesirable operational con- 
30 ditions while operating the actual machine in a region of the 
solution space defined by the Pareto frontier. 

In view of the foregoing, and with reference now to FIG. 8, 
a flow diagram 400 for use with system 300 is illustrated. In 
an off-line state 405, an exhaustive set of input parameters, 
35 such as flight envelope parameters and engine measurements 
for example, 410 are provided. At block 415, off-line multi- 
objective optimization is performed using simulator 100 and 
the aforementioned parameters as inputs. At block 420, the 
analysis results are recorded and made available for subse- 
40 quent fast lookup. In an embodiment, such availability is 
provided by storing the results in either a lookup table or in 
parametric form 425. Referring now to section 430 of flow 
diagram 400, an on-line state of process 400 is illustrated. 
Here, the on-line portion of the process begins at block 435 
45 where safety restoring parameters are retrieved via input from 
aforementioned block 425, via input from block 440, which 
provides indication of the presence of a fault from the afore- 
mentioned reasoner, via input from block 445, which pro- 
vides a current set of input parameters such as may be 
50 received from sensors 130 for example, and via input from 
block 450, where performance constraints, logistics con- 
straints, and scheduling constraints, are provided by the user. 
At block 455, the desired safety restoring parameters are 
applied at block 460 where controller 125 serves to restore 
55 actual machine 127' to safe operating conditions by using 
information from blocks 440, 445 and 450, and by applying 
multi-objective optimization analysis, as discussed above, 
from block 470. In section 465 of process 400, which illus- 
trates a portion of process 400 applicable to either on-line or 
60 off-line conditions, information from blocks 440, 445 and 
450 are provided to block 470 where multi-objective optimi- 
zation is performed in accordance with methods disclosed 
herein, and the results are provided to either block 460, as 
discussed above, or to block 475, which represents optimal 
65 results of the multi -objective optimization process. 

In an embodiment, a computer readable medium 315 for 
multi-objective HPC or HPT fault accommodation is pro- 
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vided, the computer readable medium having computer 
executable instructions for facilitating an embodiment of the 
aforementioned method. 

While an embodiment of the invention has been described 
employing an aircraft engine and aircraft engine controller, it 
will be appreciated that the scope of the invention is not so 
limited, and that the invention may also apply to any machine 
or complex machinery having a controller for controlling the 
machine or machinery. 

An embodiment of the invention may be embodied in the 
form of computer-implemented processes and apparatuses 
for practicing those processes. Embodiments of the invention 
may also be embodied in the form of a computer program 
product having computer program code containing instruc- 
tions embodied in tangible media, such as floppy diskettes, 
CD-ROMs, hard drives, USB (universal serial bus) drives, or 
any other computer readable storage medium, such as read- 
only memory (ROM), random access memory (RAM), and 
erasable-programmable read only memory (EPROM), for 
example, wherein, when the computer program code is 
loaded into and executed by a computer, the computer 
becomes an apparatus for practicing embodiments of the 
invention. Embodiments of the invention may also be embod- 
ied in the form of computer program code, for example, 
whether stored in a storage medium, loaded into and/or 
executed by a computer, or transmitted over some transmis- 
sion medium, such as over electrical wiring or cabling, 
through fiber optics, or via electromagnetic radiation, 
wherein when the computer program code is loaded into and 
executed by a computer, the computer becomes an apparatus 
for practicing embodiments of the invention. When imple- 
mented on a general-purpose microprocessor, the computer 
program code segments configure the microprocessor to cre- 
ate specific logic circuits. A technical effect of the executable 
instructions is to adjust control settings of a controller to 
accommodate for HPC and HPT faults within an aircraft 
engine. 

While the invention has been described with reference to 
exemplary embodiments, it will be understood by those 
skilled in the art that various changes may be made and 
equivalents may be substituted for elements thereof without 
departing from the scope of the invention. In addition, many 
modifications may be made to adapt a particular situation or 
material to the teachings of the invention without departing 
from the essential scope thereof. Therefore, it is intended that 
the invention not be limited to the particular embodiment 
disclosed as the best or only mode contemplated for carrying 
out this invention, but that the invention will include all 
embodiments falling within the scope of the appended claims. 
Also, in the drawings and the description, there have been 
disclosed exemplary embodiments of the invention and, 
although specific terms may have been employed, they are 
unless otherwise stated used in a generic and descriptive 
sense only and not for purposes of limitation, the scope of the 
invention therefore not being so limited. Moreover, the use of 
the terms first, second, etc. do not denote any order or impor- 
tance, but rather the terms first, second, etc. are used to dis- 
tinguish one element from another. Furthermore, the use of 
the terms a, an, etc. do not denote a limitation of quantity, but 
rather denote the presence of at least one of the referenced 
item. 
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45 What is claimed as new and desired to be protected by 
Letters Patent of the United States is: 

1. A method for multi -objective fault accommodation 
using predictive modeling, the method comprising: 

using a simulated machine that simulates a faulted actual 

50 machine; 

using a simulated controller that simulates an actual con- 
troller, the simulated machine being controlled by the 
simulated controller, and the actual machine being con- 
trolled by the actual controller; 

55 performing a multi-objective optimization process, based 

on specified control settings for the simulated controller 
and specified operational scenarios for the simulated 
machine controlled by the simulated controller, to gen- 
erate a Pareto frontier-based solution space relating per- 

60 formance of the simulated machine to settings of the 
simulated controller, including adjustment to the opera- 
tional scenarios to represent a fault condition of the 
simulated machine; and 

adjusting control settings of the actual controller, repre- 

65 sented by the simulated controller, for controlling the 
actual machine, represented by the simulated machine, 
in response to a fault condition of the actual machine, 
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based on the Pareto frontier- based solution space, to 
maximize desirable operational conditions and mini- 
mize undesirable operational conditions while operating 
the actual machine in a region of the solution space 
defined by the Pareto frontier. 5 

2. The method of claim 1, wherein: 

the actual machine is an aircraft engine. 

3. The method of claim 2, wherein: 

the fault condition is a non-catastrophic high pressure com- 
pressor (HPC) or high pressure turbine (HPT) fault con- 1 0 
dition. 

4. The method of claim 2, wherein: 

the adjusting control settings comprises optimizing actua- 
tor gains, controls modifiers, and schedules of a compo- 1 5 
nent level model of the predictive model. 

5. The method of claim 2, wherein: 

the adjusting control settings serves to reduce module tem- 
perature peaks, and to increase stall margins of the air- 
craft engine. 20 

6. The method of claim 1, wherein: 

the faulted condition of the simulated machine is represen- 
tative of the non-catastrophic faulted condition of the 
actual machine. 

7. The method of claim 1, wherein: 25 

the maximizing and the minimizing are performed off-line 

and solutions are precompiled into lookup tables that are 
stored onboard such that a solution is accessible within a 
subsecond time frame sufficient to establish safety-re- 
storing operating points for the actual machine, thereby 30 
accommodating the fault condition of the actual 
machine. 

8. The method of claim 7, wherein: 

the maximizing and the minimizing are further performed 
in a time frame greater than a subsecond subsequent to 35 
establishing the safety -restoring operating points to 
enable calculation of trade-offs on the Pareto frontier to 
incorporate contemporaneous demands on the actual 
machine and external demands. 

9. The method of claim 8, wherein: 40 

the contemporaneous demands on the actual machine 

includes one of operational safety thresholds of the 
actual machine, duration of operation of the actual 
machine, deterioration of the actual machine, risk of 
secondary damage, logistic considerations for the actual 45 
machine, planning considerations for the actual 
machine, or combinations thereof. 

10. The method of claim 1, wherein: 

the using and the adjusting are performed on-board the 
actual machine. 50 

11. The method of claim 10, wherein: 

the using mad the adjusting are performed at may time 
between consecutive operations of the actual machine. 

12. The method of claim 1, wherein: 

the adjusting control settings comprises applying a bi- 55 
directional filtering algorithm to facilitate smoothness in 
derived schedule surfaces. 

13. The method of claim 1 wherein: 

the performing a multi -objective optimization process is 
carried out in a federated manner such that subsystems 
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find their own solutions, and subsystem intercorrela- 
tions are controlled at a system level. 

14. A system for multi -objective fault accommodation 
using predictive modeling, the system comprising: 

a simulated machine that simulates a faulted actual 
machine; 

a simulated controller that simulates an actual controller, 
the simulated machine being controlled by the simulated 
controller, and the actual machine being controlled by 
the actual controller; 

a processor that performs a multi-objective optimization 
process, based on specified control settings for the simu- 
lated controller and specified operational scenarios for 
the simulated machine controlled by the simulated con- 
troller, to generate a Pareto frontier-based solution space 
relating performance of the simulated machine to set- 
tings of the simulated controller, including adjustment to 
the operational scenarios to represent a fault condition of 
the simulated machine; and 

an adjuster portion that adjusts control settings of the actual 
controller, represented by the simulated controller, for 
controlling the actual machine, represented by the simu- 
lated machine, in response to a fault condition of the 
actual machine, based on the Pareto frontier-based solu- 
tion space, to maximize desirable operational conditions 
and minimize undesirable operational conditions while 
operating the actual machine in a region of the solution 
space defined by the Pareto frontier. 

15. The system of claim 14, wherein: 

the actual machine is an aircraft engine. 

16. The system of claim 15, wherein: 

the fault condition is a non-catastrophic high pressure com- 
pressor (HPC) or high pressure turbine (HPT) fault con- 
dition. 

17. A computer readable storage medium for multi-objec- 
tive fault accommodation using predictive modeling, the 
computer readable medium comprising computer executable 
instructions which when executed on a computer perform a 
process comprising: 

simulating a faulted actual machine; 

simulating an actual controller that controls the actual 
machine; 

performing a multi-objective optimization process, based 
on specified control settings for the simulated controller 
and specified operational scenarios for the simulated 
machine controlled by the simulated controller; and 

generating a Pareto frontier-based solution space relating 
performance of the simulated machine to settings of the 
simulated controller, including adjustment to the opera- 
tional scenarios to represent a fault condition of the 
simulated machine; 

wherein the Pareto frontier-based solution space is capable 
of being used for adjusting control settings of the actual 
controller for controlling the actual machine, in 
response to a fault condition of the actual machine, to 
maximize desirable operational conditions and mini- 
mize undesirable operational conditions while operating 
the actual machine in a region of the solution space 
defined by the Pareto frontier. 





